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Abstract 



O 

' We describe a new effect which produces baryons at a first order electroweak 

m 

^ ' phase transition. It operates when there is a CP-violating field present on 

[ propagating bubble walls. The novel aspect is that it involves a purely clas- 

00 

f-*) ' sical force, which alters the motion of particles across the wall and through 

0^ ' diffusion creates a chiral asymmetry in front of the wall. We develop a tech- 

Q^' nique for computing the baryon asymmetry using the Boltzmann equation, 
and a fiuid approximation which allows us to model strong scattering effects. 



The final formula for the baryon asymmetry has a remarkably simple form. 

The last ten years have seen a growing realisation that the standard electroweak the- 
ory satisfies Sakharov's conditions for baryogenesis: B violation, departure from thermal 
equilibrium, C and CP violation [|T|. Nonperturbative B violating processes involving the 
electroweak SU{2)l chiral anomaly, appear unsupressed at high temperatures 0. The elec- 
troweak transition is weakly first order for light Higgs masses, Mh < 80 GeV 0] . It proceeds 
via bubble nucleation, with departures from thermal equilibrium on and around bubble walls 
0]. Finally, C and CP are violated, the latter via the phase in the KM matrix. The CP 
violation in the minimal theory is very small, but there may be amplification mechanisms 
which enhance it, or additional Higgs fields with CP violation in the Higgs potential (for 



1 



reviews see |5|], [|b|). 

In this Letter we study baryogenesis due to a CP violating condensate on bubble walls. 
Several mechanisms have already been pointed out through which such a condensate can 
produce a baryon asymmetry. It couples via a term in the effective action to bias the 
winding of the gauge and Higgs fields 0]. It also acts to bias hypercharge violating particle 
interactions, which in a certain constrained thermal equilibrium favors baryon production - 
'spontaneous' baryogenesis pf, 0, ||10|! IHI- 

In addition to these local effects, particle transport can carry the CP violation present 



on the wall into the unbroken phase, where the B violation rate is maximal |]12[. Until 
recently, this nonlocal barogenesis was thought to necessarily involve quantum mechanical 
interference - the idea was that the condensate causes CP violation in particle-wall scattering 
amplitudes, leading to a chiral flux (i.e. more t^'s than t^s) being injected into the unbroken 
phase where it drives baryon production. Top quarks are the obvious candidate because of 
their large mass and thus strong coupling to the wall. However significant asymmetries can 
only be produced for very thin walls because the quantum interference effect tends to be 
destroyed by the strong (QCD) scattering of the quarks, and is also WKB suppressed for 
walls much thicker than the inverse top mass. Partly for these reasons we considered tau 



leptons as an alternative because they are much more weakly coupled to the plasma [|13 
0- 

In this Letter we discuss a new, purely classical mechanism through which nonlocal 
baryogenesis can be driven. It does not rely on quantum mechanical interference, and thus 
may be calculated from a Boltzmann equation, which we shall solve analytically in a fluid 
approximation. The physical picture is extremely simple: the classical force drags an excess 
of chiral charge onto the wall, leaving a compensating deficit of chiral charge in front of 
the wall, which drives baryogenesis. Particle transport is the key to this mechanism — if 
particles are free to diffuse in the medium, they are free to respond to the chiral force on 
the wall. Conversely, if particles cannot move relative to the plasma (e.g. top quarks in the 
limit of large as) the whole effect goes away, and only local baryogenesis is possible. 



The qualitative criterion for efficient transport is that a particle should be able to diffuse 
a distance x > L, the wall width, in the time the wall takes to pass t = L/v^, with v^, the 
wall velocity. Setting ~ Dt, with D the diffusion constant, we find 

• Condition 1. < D / L (for efficient transport) 

Secondly we require that the phase space density approach a local thermal equilibrium 
(LTE) form, in which a chiral charge builds up on the wall. This requires that the equilibra- 
tion time T be smaller than the time of passage of the wall L/vw We will discuss below how 
we determine this timescale r within our calculation and find that r ~ 2)D. We therefore 
require 

• Condition 2. < L/3D (for LTE to estabhsh) 

And finally we demand that the semiclassical (WKB) approximation must be accurate. 
For this to be true, the effect should come from particles with typical momenta Pz such that 

• Condition 3. \pz\ >> L^^ (for classicality) 

Conditions 1-3 are plausibly met in the standard electroweak model, and minimal ex- 
tensions. As long as they are, we shall see that the final baryon asymmetry is to a good 
approximation independent of f^, L, and D. 

One loop calculations yield L ~ 20 — 40T^^ and in our mechanism the particles that 
dominate the effect have \pz\ ~ T, so Condition 3 is easily satisfied. We have estimated 
Tl|] D ~ 5T^^ for quarks. Calculations of the wall velocity are difficult @], but indicate 



velocities in the range ~ -05 — 1. Conditions 1 and 2 are therefore satisfied for a large 
part of the parameter space indicated by studies of the phase transition. However, we expect 
particle diffusion in front of the wall to become negligible if > Vs, the speed of sound in 
the plasma, in which case local baryogenesis should dominate. 

Recently, we pointed out in that transport effects could be important in 'spontaneous' 
baryogenesis ||^, and raised doubts about this mechanism because transport effects spoil the 
constraints imposed in local equilibrium calculations |T0[, unless the walls are very thick. 



Subsequently, one of us pointed out that transport phenomena could actually enhance this 
mechanism and this has been independently explored using a diffusion equation in [jl6 . 



We believe that our procedure provides a more complete framework within which both this 
and the classical force effect we focus on here, can be computed and give a detailed treatment 



of both effects in 11 



We begin with a derivation of the chiral force. The Lagrangian describing a fermion 
moving in the classical background of a bubble wall with CP violating condensate is 

= "^l^Kd^ - zgAZ^^')^ - mm (1) 



where QaZ^ = QaZ^^ - ^:;;iT^d^9, Qa = \\l gl + gl for t quarks The two contributions 



the CP-odd scalar field 6 which is the relative phase of the two Higgs fields in a two Iliggs 
theory, (pl'^i = Re^^, and the condensate discussed in |jl^, which may be present even in 



the minimal theory. The notation Gl implies that this is the gauge invariant combination of 
the gauge fields and Higgs phases which diagonalises the Higgs kinetic terms. We treat the 
wall as planar, and assume it has reached a stationary state in which the Higgs and gauge 
fields are functions of 2; — v^t. In this case, the field is pure gauge. 

The axial coupling in (|I|) leads to a classical chiral force as follows. In the rest frame of 
the wall Z^ = {0,0,0, Zz{z)). From the corresponding Dirac equation, setting ip oc e"*^-^, 
one finds the following dispersion relation |TT[ 



E 



S. = ±1 (2) 



P± + [yvl + gAZz 

for both particles and antiparticles. 5*^ is the component of the spin in the z direction, 
measured in the frame in which p± vanishes. In the WKB approximation, this dispersion 
relation accurately describes particles as they move across a bubble wall - the local eigen- 
states in shall form the basis of our treatment. The particles we are most interested in 
for baryogenesis are left handed particles (e.g. ^l's), and their antiparticles (ti's, which are 
right handed), because these couple to the chiral anomaly. For large \pz\i these are easily 
identifiable in terms of the eigenstates in (|^). Note that they couple oppositely to the Z 
field. 

The group velocity of a WKB wavepacket is determined from the dispersion relation by 
Vz = z = dE/dpz, and energy conservation E = = z{dE/dz) + pz{dE/dpz) then implies 



that pz = —dzE. These are of course Hamilton's equations for the motion of a particle. 
From these it is straightforward to calculate the acceleration 

where E and p± are constants of motion. 

This chiral force provided by the Zz field effectively produces a potential well which draws 
an excess chiral charge onto the wall, and leads to a compensating deficit in a 'diffusion tail' 
in front of the wall. There is net baryon production because B violation is suppressed on 
the bubble wall. 

We now seek to describe the particle excitations with dispersion relations (0) as classical 
fluids. We focus on particles with large \pz\ ^ T » m for three reasons: they dominate 
phase space, the WKB approximation is valid, and the dispersion relation simplifies so one 
can identify approximate chiral eigenstates. The Sz = +\, Pz < ^ branch, and the Sz = — |, 
Pz > branch constitute one, approximately left handed fluid L, and the other two branches 
an approximately right handed fluid R. 

The Boltzmann equation is: 

dj ^ dj + zdzf+PzdpJ = -C{f) (4) 

where z and Pz are calculated as above, and C{f) is the collision integral. This can in 
principle be solved fully. However to make it analytically tractable we truncate it with 
a fluid approximation which we now discuss. When a collision rate is large, the collision 
integral forces the distribution functions towards the local equilibrium form 

where T = v and /i are functions of z and t, and 7 = 1/(1 — f^)^/^. These parametrise 
the fluid velocity v, number density n and energy density p. We are going to treat the 
approximately left-handed excitations L and their antiparticles L as two fluids, making an 
Ansatz of the form (^) for each. 



As mentioned, the Ansatz does allow us to describe perturbations in the energy 
density, number density and velocity of each fluid, and we expect it to give a reasonable 
qualitative description of the true phase space density perturbations. As far as the temper- 
ature and velocity perturbations are concerned, we probably cannot expect this form to be 
more than qualitatively correct, because the dominant interactions which establish thermal 
equilibrium are those with the background plasma, which are also responsible for setting 
5v and 5T to zero. This is not true however of the chemical potential perturbation, which 
is only attenuated by slower chirality changing processes. So as long as we can check that 
5T/T and 5v are small, compared to yu/T, we believe that (§) should actually provide an 
accurate parametrisation of the phase space density. 

The collision integrals are evaluated in the approximation that the particle interactions 
are local, by using the Dirac spinors appropriate to the local value of m and Z^, taken to 
be constant. This is reasonable for the two body scattering effects we consider, because the 
QCD interactions are short ranged, the Debye screening length m";^^^ being smaller than 
and very much smaller than L. 

The ansatz (^) has three arbitrary functions, which are fully determined from three 

3 3 3 

independent moments of the Boltzmann equation: we take J d p, J d pE, and / d ppz- For 
a single interacting fluid these yield the continuity, energy and momentum equations. We 
are interested in the chiral density, the difference between L and L chemical potentials, since 
this quantity drives the baryon asymmetry. We work to first order in and f^. 

The fluid equations for particle minus antiparticle perturbations {6T = ST{L) — 6T{ L), 
fi = fi{L) — fi{L), 6v = 6v{L) — 6v{ L)) are, in the rest frame of the wall 

- Vy,— + -6v - aVu,—= -r^(— ) - r^(;^) (6) 

Jo <J Jo Jo Jo 

-vu,— + -5v - v^b—= -Vt— - r„(— ) - r (7) 



To 3 To ^ To '^^Tc 







— + 6— - 2cv^ -3 = -Yy^v (8) 

Jo Jo Jq 

where the shifted chemical potential difference is /i = fi — 2vgAZz, (/i) denotes the signed sum 
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of chemical potentials for particles participating in the reaction, A = (/x) = (/t — 2vwgAZz), 
is the difference between shifted L and L potentials, and prime denotes dz, a = 'k'^/27(s, 
b = noTo/po, c = In 2/14(4, C4 = vrV90, no = SCsT^^/An"^, po = 21CiT^/8'K^ and C is the 
Riemann (^-function. The derivation of these equations is simplest if one shifts the canonical 
momentum to kz = Pz + QaZz and the chemical potential to p. In this way the correct 
massless limit emerges as one expands in powers of Zz- The relevant collision integrals may 
be calculated at zero background fields flSl. is simply related to the diffusion constant 
D - it is easily seen that D = (noTo/4apo)r^^ ^ jT-K In fact we find Tt ~ |r„, ^ T/20 
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and are derived from hypercharge conserving chirality flip processes, such as those 
involving external Higgs particles. In this case, the Zz contribution to sum of chemical po- 
tentials vanishes. and F* are the rates for hypercharge violating chirality flip processes, 
which are suppressed and for these the Zz contribution does not cancel. These latter 
are the terms driving 'spontaneous' baryogenesis, in its new 'diffusion-enhanced' form. This 



enhancement was mentioned in a talk by one of us and explored independently and 
in much greater detail by Cohen et. al. |T6[. In the formalism represented by the above 
equations, the 'spontaneous baryogenesis' and 'classical force' driving terms are both in- 
cluded - this is fully discussed in ref . [ITT] . Here we focus on classical force term, which since 
the equations are linear can be considered independently of the 'spontaneous baryogenesis' 
source terms. 

We proceed to solve equations (H) to find the perturbations produced by the force term. 
We simplify by setting F^, F* , F^, and F* equal to zero. We show in [|lT| that the suppression 
they produce of the classical force effect is, for a large range of parameters, a factor between 

1 and 1. With this simplification we can derive our result in a few lines. First, from (0) we 
see that if v^jD/L < 1 the temperature fiuctuation is smaller than the velocity or chemical 
potentials by this factor (using F^ ~ 1/12D). This explains how we arrived at our Condition 

2 above. Then from (0) we find a simple relation 6v ~ f^/i. So we do indeed find the 
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temperature and velocity perturbations are small. The reason the velocity perturbation is 
small is quite general - from the continuity equation it follows that the velocity perturbation 
required to create a given chiral excess in the wall frame is proportional to the wall velocity. 
This should to be true independently of the detailed form of the phase space density, which 
as mentioned before we cannot expect to be exact. Since 5v is small, we can drop the r.h.s. 
of (H) because on the wall it is of order Vu,L/D compared to the fi term, which is small by 
Condition 1. 

We are left with a relation between the chemical potential fi and the force term on the 
wall, from fBI): 



2 In 2 QAZzm^ 

/i = — :r;—Vw — 7^ — on the wall (9) 



We can now determine (x in front of the wall as follows. Integrating (|^ and (P) (with all 
r^'s zero), gives J^^6T ^ and = 0. Then integrating (H) twice we find /^/i = 

i.e. no net integrated chemical potential perturbation is generated. This means that the 
chemical potential generated on the wall is compensated by an opposite chemical potential 
off it. As mentioned above, off the wall the equations for /i reduce to the diffusion equation, 
and it is straightforward to see that in the absence of particle number violation the only 
nontrivial solution for /i is a diffusion tail in front of the wall. This is where the chiral charge 
deficit occurs, which drives baryogenesis. 

Hence the integral of the chemical potential in front of the wall {z > 0) equals: 

/■oo 2 In 2 r 

/ dzfi = -——Vw / dzgAZ^m^ (10) 

Jo o(,3 Jwall 

Now, using the standard formula for baryon number violation 

3 r 

riB = -Vuiu'b = -^Nc^i^tL - ^J'tJ (11) 

where Tg = ^(a^T)^ is the weak sphaleron rate in the unbroken phase, is the number of 
colors, K G [0.1, 1] [0. We have re-expressed fi in terms of top quark and antiquark chemical 
potentials. We arrive at a formula for the baryon to entropy ratio: 

8 



ub 1351n2Ka^ f ^jtl^qaZ^ ^^na^, f ^jnlgAZ^ 



, dz ' ^ 4 J / dz ' (12) 

where s = '^g^T^ is the entropy density, with g^ ^ 100 the effective number of degrees of 
freedom. 

This result is remarkably simple - all dependence on the wall velocity, thickness and 
the diffusion constant drops out, provided Conditions 1 and 2 are satisfied. It is also quite 
large: (mt/T)^ ~ 1, so ns/s ~ 4 x IO^^kOcp where 6cp characterises the strength of the 



CP violation. In a longer paper |Tl| we give a more detailed derivation of ([T2|) with a full 
discussion of parameter dependences, including the effect of the terms we have neglected 
here. 

The dependence in ( |T2D means that, at least with standard model-like Yukawa cou- 
plings, the top quark dominates the effect. The mass-over-temperature suppression can be 
significant, if the Zz field is localised on the front of the wall where the Higgs vev is small. 

The calculation of the classical force effect above uses the opposite (WKB) approximation 



to those employed in quantum mechanical reflection calculations (thin walls) [0, |13]). The 
classical force calculation is in some respects 'cleaner', because the production of chiral charge 
and its diffusion are treated together. The classical force affects particles from all parts of 
the spectrum, mostly with typical energies E ^ T, and with no preferential direction, while 
the quantum mechanical effect comes mainly from particles with a very definite ingoing 
momentum perpendicular to the wall: Pz ~ rriH (Higgs mass). The quantum result falls off 
strongly with L (at least as L^"^) as the WKB approximation becomes good. The quantum 
result also has a dependence coming from the diffusion time in the medium, which the 
classical result loses because the force term is proportional to Vw 

Finally, we mention possible extensions of these methods. In the above treatment we 
have comletely ignored collective plasma effects - Debye screening. Landau damping etc, 
and merely treated local particle interactions. The Boltzmann equation is easily modified to 
include these effects, with force terms due to the electric (and magnetic) fields, which are 
solved for self-consistently. The fluid truncation may be a useful way to compute the bubble 



wall velocity (at least the friction due to top quarks), and we shall return to this in future 
work. We intend also to extend these methods to study Z-condensation in the standard 
model in the presence of strong interactions. 
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